Length scales coupling for nonlinear dynamical problems in magnetism 
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The dynamics of real magnets is often governed by several interacting processes taking place si- 
multaneously at different length scales. For dynamical simulations the relevant length scales should 
be coupled, and the energy transfer accurately described. We show that in this case the micromag- 
netic theory is not always reliable. We present a coarse-graining approach applicable to nonlinear 
problems, which provides a unified description of all relevant length scales, allowing a smooth, seam- 
less coupling. The simulations performed on model systems show that the coarse-graining approach 
achieves nearly the precision of all-atom simulations. 

PACS numbers: 75.40.Mg, 75.60.Ej, 75.30.Hx, 02.70.Dh 



The dynamics of magnetization in real magnets is often 
governed by defects: impurities, vacancies, grain bound- 
aries, etc. They serve as nucleation centers for magne- 
tization reversal or as pinning centers for domain walls 
(DW), due to very localized (several lattice constants) 
but significant changes in the magnetic interactions. As 
a result, the torques acting on individual atomic spins 
in the vicinity of the defect are strongly non-uniform 
at the atomic scale. These small-scale inhomogenieties 
can result, however, in large-scale (hundreds or thou- 
sands of angstroms) changes of the magnetization dis- 
tribution (e.g. nucleation of a domain). At larger length 
scales (microns) , the magnons propagate into the bulk of 
the sample carrying away significant energy. Such mul- 
tiscale phenomena are encountered in the depinning of 
DW pL ||, ||, Q, the influence of surfaces and interfaces 
on the spins in the bulk @, [| fj], ||] , etc. These phenomena 
determine basic magnetic properties such as the coercive 
field and the dynamics of magnetization reversal. 

Accurate simulations of such phenomena might be pos- 
sible if all the interactions and dynamics could be treated 
at the atomic scale Q , but for a sample of the size of few 
tens of microns, the number of spins becomes much too 
large for modern supercomputers. Atomistic modeling 
can be used, however, for the regions close to the defects. 
For the large-scale processes occuring far from the defects 
a collective description should be used, e.g. micromag- 
netic (MM) theory. To simulate accurately the multiscale 
processes, all the relevant scales should be modeled simul- 
taneously, and seamlessly coupled to describe correctly 
the energy transfer. The simplest solution would be to 
use the standard MM approach with gradually increasing 
resolution (reducing the size of the computational cell) 
near the defect, until every micromagnetic cell contains 
a single atomic spin. At this ultimate resolution, the mi- 
cromagnetic Landau-Lifshitz-Gilbert equation coincides 
with the equation of precession of an individual atomic 
spin. 

In this paper, using simulations performed on model 
1-D chain systems, we show that such an approach is 



often inadequate for length scale coupling, giving results 
which differ drastically from the exact solutions obtained 
by all-atom simulations. Thus, for an important class of 
dynamical problems in magnetism, a new approach is 
needed. Here, we present a theoretical framework along 
with a computational scheme, which allows treatment of 
all length scales in exactly the same manner, enabling a 
truly seamless linking of different scales. The proposed 
approach is applicable to essentially nonlinear problems, 
and yields reliable results which are in excellent agree- 
ment with exact atomic simulations. 

Length scale coupling has been extensively discussed in 
the context of lattice dynamics modeling O, |l5f| 
and successfully applied for the study of crack prop- 
agation and nanoindentation |b|, Analo gous 
schemes can be constructed for linear magnetic prob- 
lems |fr7| , but most of the interesting problems in mag- 
netism are essentially nonlinear due to the constraint 
Si + Sy + S 2 = S 2 (S x . VtZ are the components of the 
spin and S is its length; everywhere below we assume 
5=1). Moreover, the lattice constant a is the relevant 
length scale for a typical problem of lattice dynamics, 
while magnetic problems have one more intrinsic length 
scale, the DW width A. Therefore, magnetic multiscale 
simulations present some unique problems. 

Our approach to linking length scales in dynamical 
magnetic simulations is based on statistical coarse grain- 
ing (CG), which has been successfully used before |l2|, [l7|] 
to define suitable collective variables ("gross variables"). 
For simplicity, we consider only the case of zero tem- 
perature, thus restricting the analysis to purely dynam- 
ical problems, with dissipation and thermal noise effects 
absent. We assume that the system under considera- 
tion is a ferromagnet made of identical classical spins S^, 
(|S M | = 1), located at the /i-th site of the crystalline lat- 
tice (greek indices enumerate the atomic lattice sites). 
The system is described by a rather general spin Hamil- 
tonian 

h = h° + v, n° = J2K = J2 J ^ S » S »> w 
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where Ti° describes the isotropic exchange interaction, 
while V <C Ti° represents all the other, much weaker 
interactions, and J M „ is the exchange interaction be- 
tween the atomic spins at the sites fj, and v. We will 
need a Hamiltonian formulation of the magnetization 
dynamics, therefore below we describe spin by conju- 
gate canonical variables ai M = 2 sin (0 M /2) cos M and 
a2fi = 2 sin (6*^/2) sin M , where 9 a and <fi a are the az- 
imuthal and the polar angles of the spin vector corre- 
spondingly. We assume that for the simulations of the 
large-scale regions located far from the defects, the finite- 
element method (FEM) is used. Thus, we define the 
gross variables, which describe the state of the system 
at a given length scale, by averaging the atomic degrees 
of freedom: ay = ^ U,j and a ?3 = J2u U,j a 2 M , 
where the index j corresponds already to the computa- 
tional FEM nodes. The weight function localized 
near the node j must satisfy the normalization condition 
fuj = 1. The choice of this function was discussed 
for linear magnetic multiscale problems; our expe- 
rience shows that the following piecewise-linear choice is 
sufficient: f^j = /o|/i-/ij|/l/i,' _ ^j-i| for fi e [fij-x, 
where \ij is the atomic position of the j-th computational 
node and fo is the normalization. 

The main assumption of the approaches based on the 
coarse graining ideology (such as the non-equilibrium sta- 
tistical operator Q or projection operator Q methods) 
is that the atomic degrees of freedom are in local equi- 
librium, and the equilibrium conditions are determined 
by the gross variables ^o). This means that all atomic 
spins located near the node j are nearly parallel to the 
direction defined by the the variables cti,2ji the exchange 
part Ho of the Hamiltonian ([!]) is responsible for that 
(thus, the coarse-graining approach can be used also for 
finite temperature if it is small in comparison with the 
Curie temperature). Then, the distribution function for 
the atomic variables is 

p = Q^expi-pH'), (2) 

where Q is the statistical integral, and the torques Fj 
and Gj are determined from the conditions of local equi- 
librium: Yipfnj ("im) = and ( a 2 M ) = a> 2j , 
where the angular brackets mean the statistical averaging 
with the distribution function @. 

However, the exchange Hamiltonian H.q is rotationally 
invariant, and a zero- frequency mode is present in the 
statistical integral Q in Eq. |^. In the linear approxi- 
mation, this leads to a divergence, and the account of 
nonlinear effects becomes possible. Such a situation has 
been studied e.g. by Bogoliubov in the work on the small 
polaron problem, see e.g. Ref. [|lj. The general ideas of 
that work can be applied here. 

We need to separate the (nonlinear) motion of gross 



variables from the (quasilinear) motion of the atomic- 
scale variables to exclude the secular terms causing di- 
vergence of Q. To do this, we note that the charac- 
terise scale of nonlinear magnetic structures, the DW 
width A, is much bigger than the interatomic distance 
a, and the separation of the slow and the fast variables 
can be done in the spirit of the modified adiabatic the- 
ory pl| . Instead of using the parameters a\ t 2j, we de- 
fine a local coordinate frame, rotated with respect to the 
"laboratory" frame by the angles 9j and cf>j such that 
ctij = 2 sin (6j/2) cos cj>j and ct2j = 2 sin (Oj/2) sin (j>j. 
To exclude the secular terms, we use the condition 

U,j = E„ U,j = 0, where a[% are de- 

fined already in the local coordinate frame associated 
with j-th node. These conditions are local, since the 
weighting function f a j is local, and the statistical inte- 
gral Q is also calculated in the local coordinate frame. As 
a result, the equations of motion for the gross variables 
9j and <f>j can be obtained in closed form: 

jU kn 

sm6 j fa = ^/^(dV/da^-J^MjkDkna^, (3) 

M kn 

where M jk = J2ufvdfn,k, and D jk = D° jk ~ 
djSk/Sa - dkSj/So + d SjS k /SQ, where D° k = 

\J2^u fii,j^uU,kj , J&v is the Hessian matrix of the 
exchange Hamiltonian, S k is the eigenvector of D® k cor- 
responding to the zero eigenvalue, and dj = Yli^jii 
do = J2i di, So — Yli Si- Note that the matrices M and 
D should be calculated only once and used later during 
simulations. 

Now, let us focus on the numerical results obtained us- 
ing the procedure described above. The first case tests 
the ability of our CG scheme to handle the dynamics of 
an essentially nonlinear object, e.g. a domain wall. We 
consider a I-D chain of N — 465 spins with ferromagnetic 
nearest-neighbor coupling J — 25, and antiferromag- 
netic next-nearest-neighbor coupling J' = — 7J = —2.5. 
The spins possess single-ion anisotropy of easy-axis type, 
K = 0.01 (the easy axis is directed along the y-direction) ; 
the dipole-dipole interactions are neglected for simplicity. 
The ends of the chain are different from the bulk; this 
represents e.g. a nanowire which has the growth defects 
at the ends causing the spin Hamiltonian parameters J, 
J' and K to be different from their bulk values. Thus, six 
spins at one end of the chain are chosen to have the fol- 
lowing parameters: Joi = 0.5J, J12 = 0.6J, J23 = 0.7J, 
J34 = 0.8J, J 45 = J 5 6 = J,K = -OAK, K x = -0.2K, 
K 2 = 0, K 3 = 0.2K, K4, = OAK, K 5 = 0.8K. At the 
other end of the chain, the last six spins have the same 
parameters in reverse order. Initially, the chain contains 
a Neel DW, with spins rotating in the x—y plane from 
the direction along the y-axis to the opposite direction. 




FIG. 3: Time dependence of the total magnetization Stot(i) 
of the chain calculated by three methods: large white circles 
— atomic simulations, small black circles — MM simulations, 
crosses — CG scheme. Initial conditions correspond to chain 
uniformly magnetized along the j/-axis. Since the CG results 
almost coincide with the atomistic ones, the crosses appear 
mostly in the centers of the white circles. 

the absence of the defects at the ends of the chain, all 
the spins would rotate in unison. However, in the pres- 
ence of the perturbation caused by the defects, different 
spins rotate with slightly different rates, and the system's 
motion becomes stochastic. Therefore, after some time, 
the Zeeman energy is transferred to the energy of short- 
wavelength magnons, leading to a gradual decrease of the 
system's total magnetization |^3|, |24| . 

The chain's total magnetization as functions of time is 
shown in Fig. ||, with the time unit r = 0Aw/\H sin0| w 
66. Magnetization profiles at t = 60r are presented in 
Fig. |]. Again, one can see that the MM approach docs 
not describe the behavior of the system correctly even 
at a qualitative level. The CG scheme performs much 
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FIG. 4: The magnetization profile in the chain (the angles co^ 
and ipn) at t = 60r calculated by three methods: large white 
circles — atomic simulations, small black circles — MM sim- 
ulations, crosses — CG scheme. Initial conditions correspond 
to chain uniformly magnetized along the y-axis. 

better, giving the results very close to the exact atomistic 
solution. 

Such a drastic difference in the two schemes, CG and 
MM, is caused by the essential nonlinearity of magne- 
tization dynamics. Due to the presence of atomic-scale 
inhomogenieties, in the course of the system's motion, a 
considerable number of short-wavelength excitations are 
generated near the defect. Wavelengths of these magnons 
are smaller than the size of the computational cell, and 
the micromagnetic theory neglects them completely, al- 
though the collective effect of a large number of such 
modes is significant. In contrast, the coarse-graining ap- 
proach does take the atomic-scale modes into account, 
and describes very well the collective properties of an 
ensemble of these short-wavelength excitations. 

In summary, we have shown that the standard micro- 
magnetic theory does not handle correctly the dynam- 
ics of nonlinear multiscale magnetic processes. We have 
suggested another approach, based on statistical coarse 
graining, which is applicable to nonlinear problems. Nu- 
merical tests on 1-D systems show that the CG scheme, 
in contrast with the MM approach, gives almost exact 
results for rather large time spans. It is worthwhile to 
note that the only essential ingredient of the coarse grain- 
ing approach is the standard theory of local equilibrium. 
Therefore, the CG scheme can be applied to a large set 
of 2-D and 3-D problems, including dissipation, thermal 
noise, interaction with other degrees of freedom (spin- 
phonon interactions) etc. 
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